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We study the scaling properties of polymers in a d-dimensional medium with quenched defects 
that have power law correlations ~ r~ a for large separations r. This type of disorder is known to be 
relevant for magnetic phase transitions. We find strong evidence that this is true also for the polymer 
case. Applying the field-theoretical renormalization group approach we perform calculations both 
in a double expansion in e = 4 — d and 8 = 4 — a up to the 1-loop order and secondly in a 
fixed dimension (d = 3) approach up to the 2-loop approximation for different fixed values of 
the correlation parameter, 2 < a < 3. In the latter case the numerical results need appropriate 
resummation. We find that the asymptotic behavior of self-avoiding walks in three dimensions and 
long-range-correlated disorder is governed by a set of separate exponents. In particular, we give 
estimates for the v and 7 exponents as well as for the correction-to-scaling exponent 10. The latter 
exponent is also calculated for the general m-vector model with m — 1, 2, 3. 

O 

O I I. INTRODUCTION 

The influence of structural disorder on the critical behavior of various kinds of condensed matter remains one of the 

^ ■ central problems in physics. In this paper, we are interested in the scaling laws that govern the behavior of polymers 

in disordered media when the defects are correlated or belong to some porous or sponge-like structure. Our main 

' ■ question of interest will be: does a small amount of correlated quenched structural defects in the medium induce 

l» changes to the universal properties of a polymer macromolecule? 

It is well established that the universal scaling properties of long flexible polymer chains in a good solvent are 

__ ' perfectly described within a model of self-avoiding walks (SAWs) on a regular lattice B. The limit of SAWs with 

an infinite number of steps may be mapped to a formal in —> limit of the m-vector model at its critical point [0 . 

^^ ■ In particular, for the average square end-to-end distance R e and the number of configurations Z^ of a SAW with JV 

steps on a regular lattice one finds in the asymptotic limit N — ► 00: 



(R%) ~ N 2u , Z N ~ z^A^- 1 (1) 

T3 ; 

where v and 7 are the universal correlation length and susceptibility exponents for the m = model that only depend 
ji ' on the space dimensionality d and z is a non-universal fugacity. For d — 3 the exponents read [[| v — 0.5882 ± 0.0011 
and 7 = 1.1596 ± 0.0020; whereas for d — 2 exact values v = 3/4 and 7 = 43/32 are known B. 

The problem of SAWs on randomly diluted lattices, which may serve as a model of linear polymers in a porous 
medium, has been the subject of intensive discussion M @, u\ H, pi HG, prl n3, O) . A recent review on SAW statistics on 



random lattices is given in ref. 1 10 . The numerical results for these systems available from Monte-Carlo simulations, 
exact enumeration and analytical treatment also cover the non-universal properties. Nonetheless, even apart from 
the numerical values of scaling exponents the question if a given form of disorder affects the scaling behavior has not 
been settled in general. 

A frequently studied type of random lattice is the lattice that is diluted to the percolation threshold p0[. Here, 
one is interested in the behavior of a SAW on the percolation cluster. Scaling laws (|l|) hold with exponents that differ 
from their counter parts on a regular lattice at d = 2 and d — 3 P, pi ■ Apparently, this results from the fact that the 
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percolation cluster itself is characterized by a fractal dimension which differs from d, the Euclidean one. Moreover, the 
scaling of the averaged moments of Z^ @) on the backbone of a percolation cluster possesses multifractal behavior 
J12| . In our study, however, we address another type of disorder, when the lattice is well above the percolation 
threshold. In this case the dimension of the support does not change and it is not clear a priori whether the SAW 
asymptotic exponents (0) will be influenced. 

We approach this question for the case of long-range-correlated disorder using the connection between the scaling 
properties of polymers and magnets. So let us first turn our attention to the magnetic problem. While it is intuitively 
clear that strong disorder destroys the magnetic ordering, a much more subtle question is what happens at weak 
dilution by a non-magnetic component, i.e. well above the percolation threshold (weak disorder) [H. It has been 
argued Q that the presence of a non-correlated (or short-range-correlated) quenched disorder has a nontrivial effect 
on the critical behavior of magnetic systems, only if the specific heat critical exponent a of the pure magnet is positive. 
This statement is often called the Harris criterion. However, one should be careful in applying this "naively" to the 
SAW problem. Indeed, although the critical exponent a of a SAW on the d = 3 dimensional pure lattice is positive 
H (a(d = 3) = 0.235 ± 0.003), a weak quenched short-range-correlated disorder does not alter the SAW critical 
exponents. This statement has been proven by Harris |7J] and confirmed later by renormalization group results |6[. 

Note, that in the works mentioned above only uncorrelated quenched defects were investigated. In 1980-ies the 
model of a disordered d-dimensional system with so-called "extended" structural defects [If], [17| was developed. These 
defects are considered as quenched and correlated in a subspace of Ed dimensions, and randomly distributed in the 
remaining d — £d dimensions. This model may be applied for small densities of defects. The integer values of Ed have a 
direct physical interpretation: £d = corresponds to short-range-correlated point- like defects, and the cases Ed = 1, 2 
are related, respectively, to lines and planes of impurities. To give an interpretation to non-integer values of Ed, one 
may consider patterns of extended defects like aggregation clusters, and treat Ed as the fractal dimension of these 
clusters JL8[ . In this interpretation the defect patterns are fractal, while the support of the system is the complement 
of this fractal and will in general not be fractal itself. In ref. u& the critical behavior of 0(m) symmetric magnets with 
extended defects with parallel orientation was investigated evaluating the renormalization group equations by a double 
e = 4 — d, Ed expansion. It was found that the scaling is affected by these kinds of defects and the critical exponents 
were calculated in this scheme. The static and dynamic critical properties of m-component cubic-anisotropic systems 
with extended-defects were studied in the refs. Jl9(] using a double e, e' = £ + e^-expansion again finding a change of 
critical behavior when Ed is increased. 

In further work |2C[| attention concentrated on disordered systems with "random-temperature" disorder, arising from 
a small density of impurities that cause random variations in the local transition temperature T c (x). The fluctuations 
in T c (x) are characterized by a correlation function, that falls off according to a power law: ~ x~ a at large distances 
x. It was shown, that in the presence of long-range-correlated disorder the Harris criterion is modified: for a < d the 
disorder is relevant, if the correlation length critical exponent of the pure system obeys v < 2/a. An m-vector model 
of this type was evaluated using a renormalization group expansion in the parameters £ — A — d 1 <5 = 4 — a up to the 
linear approximation. An additional renormalization group fixed point corresponding to the long-range correlated 
disorder was found. In the following we will denote this as the 'LR' fixed point. The correlation-length exponent 
was evaluated in this linear approximation as v = 2/a and it was argued, that this scaling relation is exact and also 
holds in higher order approximation. However, this result was questioned recently in refs. |21| , where the static 
and dynamic properties of 3d systems with long-range-correlated disorder were studied in a renormalization group 
approach using a 2-loop approximation. There is an essential discrepancy between the latter results and those found 
from the e, ^-expansion. Nevertheless it was qualitatively confirmed by both approaches that long-range-correlated 
disorder leads to a new universality class for these magnetic systems. Note, that the variable a is a global parameter: 
together with the space dimension d and the number of components m of the order parameter it fixes the universal 
values of the critical exponents. 

While the influence of long-range-correlated disorder on the magnetic phase transition has been the subject of 
considerable interest, the effect of long-range-correlated disorder on the scaling properties of polymers remains unclear 
and is generally not considered as settled. Here, we address the question of the asymptotic behavior of polymers in 
long-range-correlated disorder with algebraically decaying correlations |2Q] . While the linear approximation of the 
double e, 5-expansion indicates qualitatively the existence of the LR fixed point for polymers, it leads to unphysical 
quantitative results |13|] . For this reason we present here an analysis of the 2-loop approximation using the fixed 
a, d-technique that leads to physically meaningful results for the scaling behavior of polymers in the LR regime. 

Our paper is organized as follows: in the following section Owe present the model, in section [II the renormalization 



procedure is discussed and we reproduce the results of the e, 5-expansion. In section IV we apply resummation 
techniques to analyse the renormalization group functions in the two-loop approximation and find that the asymptotic 
behavior of self-avoiding walks in a three-dimensional medium with long-range-correlated disorder is governed by 
a new set of exponents. For the exponents we present quantitative estimates. Section M concludes our study. 
Some additional information about the properties of magnetic phase transitions in systems with long-range-correlated 



quenched disorder is presented in the appendix. 

II. THE MODEL 

To study the universal properties of polymers in porous media with long-range-correlated quenched structural 
defects, we turn our attention to the investigation of the appropriate m-vector model in the polymer limit. We 
consider the model of an m-vector magnet, that is described by the following Hamiltonian: 



H = d d x 



\(fjM> + 8lM>(*))P + (V<?) 2 )+^) 2 



(2) 



here, 4> is an m-component field: <p — {(f) 1 ■ ■ -f/)" 1 }, /io and u o are the bare mass and the coupling of the undiluted 
magnetic model, 6fj,o(x) represents the quenched random-temperature disorder, with: 

«tyo(aO» = 0, 
-({5/j,o(x)5no(y))) = g(\x - y\), 

where ((• • •)) denotes an average over spatially homogeneous and isotropic quenched disorder. The form of the pair 
correlation function g(r) is chosen to fall off with distance according to a power law [ p0[ : 

g(r) ~ r~ a (3) 

for large r, where a is a constant. 

We consider quenched disorder and average the free energy over different configurations of the disorder. To this end 
we apply the replica method and construct an effective Hamiltonian for the m-vector model with long-range-correlated 
disorder p0[ : 

H cS = J2 / dd 4^( M0 ^ + ( ^ )2 ) + ¥ (< ^ )2 ]" S fd d xd d yg(\x-y\)$l(x)^(y). (4) 

Here, the replica interaction vertex g(r) is the correlation function given in Eq. (pto, Greek indices denote replicas and 
the replica limit n — > is implied. 

For small k the Fourier-transform g(k) of g(r) reads: 

g(k)~v + w \k\ a - d . (5) 

Note, that in the case of random uncorrelated point-like defects the site-occupation correlation function formally 
reads: g(\x — y\) ~ 5(\x — y\), and its Fourier transform obeys: 

g(k) ~ uo- (6) 

Comparing Eqs. (pj) and rtq), it is obvious that the case g{r) ~ r~ d corresponds to random uncorrelated point-like 
disorder. Moreover, different integer values of a correspond to uncorrelated extended impurities of random orientations. 
So, the correlation function in Eq. (H) with a = d — 1 describes straight lines of impurities of random orientation 
whereas random planes of impurities correspond p6| to a = d — 2. In terms of the fractal interpretation given in the 
introduction, the general case corresponds to SAWs on the complement of a fractal with dimension e^ = d — a. 

Writing Eq. (^) in momentum space and taking Eq. (g) into account, one obtains an effective Hamiltonian with 
three bare couplings uq, Vo, wq. For a > d the wo-term is irrelevant in the renormalization group sense and one obtains 
the effective Hamiltonian of the quenched diluted (uncorrelated) m- vector model J22J with two couplings uo,vq. For 
a < d we have, in addition to the momentum-independent couplings, the momentum dependent one wok a ~ d . Note 
that g(k) must be positive being the Fourier image of the correlation function. This implies wq > for small k. Also 
the coupling uq must be positive, otherwise the pure system would undergo a first order phase transition. 

The critical behavior of the model in Eq. (H) with m > 1 has been investigated (20| gl], ^3) using the renormalization 
group approach p4j . We are interested in the polymer limit m — ► of this model interpreting it as a model for polymers 
in a disordered medium. Note, that this limit is not trivial. For the case uq ^ 0, t>o 7^ 0, wq — the "naive" RG 
analysis leads to controversial results about the absence of a stable fixed point and thus to the absence of the second 



order phase transition [g). As noticed by Kim H, once the limit m, n — ► has been taken, the ito and ^o terms are of 
the same symmetry, and an effective Hamiltonian with one coupling of 0(mn = 0) symmetry results. This leads to 
the conclusion that weak quenched uncorrelated disorder is irrelevant for polymers as long as vq < uq- 

Our present analysis takes these symmetry properties into account. In the case of the Hamiltonian with a term for 
long-range-correlated disorder Eq. (0) we pass to an effective Hamiltonian p3[ with only two couplings Uq = uq — vq 
and wo (in what follows below we will keep the notation uq for this new coupling Uq). In discrete momentum space 
this effective Hamiltonian reads: 



n _. n 

n «* = E E 2 (/io + fc2)( ^ )2 + jiT, Y, s ^ + fc 2 + fc 3 + h) ($ a kl $t) ($tj a ki ) - (7) 

n 

^fE E \k\ a - d 5(k l+ k 2 + k)5{k, + k i -k)[$tJt){^ k M i ), m,n^Q. 



a/3 kk\k^k^,k^ 

Here, the <$(fc) represent products of Kronecker symbols and the notation (</></>) implies a scalar product. Note that 

the wo-term itroduces interactions between the replicas and contains the power of an internal momentum. Again, it 
may be shown that for a — d in the limit m, n — > the uq and wq terms are of the same symmetry and one is left 
with an 0(mn — 0)-vector model with only one coupling (uo — wq)- 

III. THE RENORMALIZATION 

In order to extract the critical behavior of the model, we use the field-theoretical renormalization group (RG) 
method. We choose the massive field theory scheme with renormalization at non-zero mass and zero external momenta 



25[ that leads to Callan-Symanzik ec 



(M) 



uations for the renormalized one-particle irreducible vertex functions T R . In our 
case the renormalization conditions 21j] are written both in fixed d and a. The renormalized mass /i and renormalized 
couplings u, w are defined by: 

fi 2 = r^ ) (fc,/i 2 ,u,w)| fc=0 , 
fj, 4 - d u = T)i iU ({k},iJ?,u,w)\k=o, 
^- a w = r^({fc} jAi 2 , U , W )| fc= o. 

Here, T Ru and T Rw are the contributions to the four-point vertex function T R that correspond to u- and w-term 
symmetry, respectively. Asymptotically close to the critical point the M-point renormalized vertex functions obey 
the homogeneous Callan-Symanzik equation J24J: 

|^ + E^(i^})^:-Y^(^>)} r ^ ) (W^ 2 '^}) = ' ( g ) 

here, Vi = u,w. The change of the couplings u,w under renormalization defines a flow in parametric space that is 
governed by the corresponding j3- functions fl u (u,w), f3 w (u,w). The fixed points u*,w* of this flow are given by the 
solutions of the system of equations: j3 u (u* , w*) = 0, f3 w (u*,w*) = 0. The stable fixed point is defined as the fixed 
point where the stability matrix: 

B ^ av (9) 

possesses eigenvalues u% with positive real parts. The accessible stable fixed point corresponds to the critical point 
of the system. The fixed point is accessible if it can be reached along flow lines starting from allowed initial values 
uq^Wq > 0. At the fixed point we define the correlation length and pair correlation function critical exponents v and 
f] by 

v- 1 = 2-7 (u*,O-70=(u*,O, (10) 

V = l${u*,w*), (11) 

where 7^2 is the exponent that corresponds to the two-point vertex function M = 2 with a <\? insertion. Other critical 
exponents may be obtained from familiar scaling laws. For example, for the susceptibility exponent 7 one has: 

7 = 1/(2-77). (12) 



According to the RG prescriptions given above, the RG functions are obtained in form of a series in the renormalized 
couplings. In the one- loop approximation the result reads [pL3|| : 



Pu 



u 2 h 



52uw 



1 



;^4 



+ (2S-e)-w 2 I 3 , 



(13) 



w = -6 



w 2 I 2 



- e- [wuh - w 2 h] , 



(14) 



U T X W T 



S^-I, 



70 2 =e-i\-0 — i 2 , 70 
Here, U are one-loop integrals that depend on the space dimension d and the parameter a: 



h = 



dq 



(q 2 + I) 2 



h = 



dqq a - d 
(q 2 + I) 2 



dq, q 



2(a-d) 



(q 2 



1) 



h = 



d 

5P 



dqq a 



-d 



[q + kf + l 



(15) 



(16) 



fe 2 =0 



Note that contrary to the usual 4 theory the 7^ function in Eq. ( |15| ) is nonzero already in one-loop order. This is 
due to the fc-dependence of the integral I4 in Eq. (jig). 

There are two ways to proceed in order to obtain the qualitative characteristics of the critical behavior of the model. 
One can consider the polynomials in Eqs. (|l3|), ( |14| ) for fixed a, d and look for the solution of the fixed point equations. 
It is easy to check that these one-loop equations do not have any stable accessible fixed points for d < 4. The other 
scheme to evaluate these equations is a double expansion in e = 4 — d and 6 — 4 — a as proposed by Weinrib and 
Halperin pOfl . Formerly [ |l3| , we exploited this up to the linear approximation. For completeness, we here note those 
results. Substituting the loop integrals in Eqs. (J13|)-(|l5|) by their expansion in e = 4 — d and S = 4 — a , one obtains 
the 3 fixed points given in thea Table H. We may draw the following conclusions from these first order results: Three 
distinct accessible fixed points are found to be stable in different regions of the a, d-plane. The Gaussian (G) fixed 
point, the pure (P) SAW fixed point and the long-range (LR) disorder SAW fixed point. The corresponding regions 
in the a, d-plane are marked by I, II and III in fig. QJ. In the region IV no stable fixed point is accessible. 

For the correlation length critical exponent of the SAW, one finds distinct values ^p U ro for the pure fixed point and 
i/lr for the long-range fixed point. Taking into account that the accessible values of the couplings are u > 0, w > 0, 
one finds that the long-range stable fixed point is accessible only for <5<e<2<5, ord<a<2 + d/2, a region where 
power counting in Eq. (fn) shows that the disorder is irrelevant. In this sense the region III for the stability of the LR 
fixed point is unphysical. Formally, the first order results for d < 4 read: 



^purc 



= 1/2- 
= 1/2 + 



-e/16, 
6/8, 



6 < e/2, 

e/2 < S < e. 



(17) 



Thus, in this linear approximation the asymptotic behavior of polymers is governed by a distinct exponent i/lr in the 
region III of the parameter plane a, d . 

Something similar happens if the e,<5-expansion is applied to study models of m-vector magnets with long-range- 
correlated quenched disorder. For comparison, using the results of Weinrib and Halperin |2(| we get the phase diagram 
presented for m = 2 in the fig. 0. Although the critical behavior of the long-range-correlated universality class appears 
there for a < d where it is relevant by power counting (region III in the fig. m) this region in the d— a plane is separated 
from the critical behavior of the 0(m = 2) universality class by a region, where no accessible stable fixed points are 
present (IV in fig. pi). It means that also in the case of magnets, as well as for polymers the first order e, ^-expansion 
leads to a controversial phase diagram (compare figs. El and ph. So our first order results should be considered as 
purely qualitative and in order to obtain a clear picture and more reliable information, we proceed to higher order 
calculations. 



IV. THE RESUMMATION AND THE RESULTS 



Fortunately, to investigate the 2-loop approximation in a fixed d and a approach we need not recalculate the 
intermediate expressions of perturbation theory for the vertex functions. Instead, we may make use of m — > limit 
of the appropriate m-vector model, investigated recently |2l| . Starting from the two-loop expressions of ref . |21[ for 
the RG functions of the m-vector magnet with long-range-correlated disorder and applying the symmetry arguments 



II [13| for the polymer limit m = as explained in Section || we get the following expressions for the d — 3 RG 
functions of the model in Eq. (0): 

95 1 

u (u,w) = -u + u 2 - (3/i(o) - f 2 {a))uw- — u 3 + -b 2 (a)u 2 w - 

(63(a) + ~& 6 (a))uu; 2 + / 3 (a)w 2 + 6 5 (a)^ 3 , (18) 

7/7/1 x, rs 

/9u,(u, u>) = -(4 - a)w - (/i(a) - f 2 (a))w 2 + — + &i (a)w 3 - — — u u> + -6i 2 (a)ww 2 (19) 

2 21o 4 

7^ (it, «;) = -f 2 (a)w + YQg" 2 + ci(a)w 2 - -c 2 (a)uw, (20) 

702 (m,w) = -m- -fi(a)w- —u 2 - c 3 (a)w 2 + -c 4 (a)uw. (21) 



Here, the coefficients /j(a) are expressed in terms of the one- loop integrals in Eq. (16), 6j(a) and Ci{a) originate 



from twodoop integrals and are tabulated in ref. pi| for d = 3 and different values of the parameter a in the range 

3m r ,,, . 3w 



2 < a < 3. The series are normalized by a standard change of variables u — » ^Ji, iu — ♦ ^/i, so that the coefficients 



in (3 U become 1 in modulus. 

The RG functions listed above have the form of a divergent series, with zero radius of convergence p6[ , familiar to 
the theory of critical phenomena p4fl . If the nature of the divergence is such that the series are asymptotic, then the 
situation is, at least in principle, controllable: in this case a good estimate for the sum of the series is obtained by 
keeping a certain number of the first terms ("optimal truncation") or applying an appropriate resummation procedure. 

For the case of the pure 3-dimensional 4 theory it is known that the perturbation series are asymptotic, and Borel- 
summability in 3 dimensions has been proven p7| . The situation of the random-site Ising model is less satisfactory 
than for the pure system [14fl . For instance, the asymptotic parameter in the disordered system is y/e instead of e, 



and the /3-functions, computed at two loops show no stable fixed points. Bray et al |28 and McKane [g9[ studied the 
asymptotic expansion for the free energy of the random-site Ising model in the zero-dimensional case, and the model 
was found to be non-Borel summable. However, recently pCfl the Borel summability of the perturbation expansion 
for the zero-dimensional disordered Ising model was proven analytically, provided that the summation is carried out 
in two steps: first, in the coupling of the pure Ising model and subsequently in the variance of the quenched disorder. 
In our case, the summability of the series in Eq. (118) is open. Nevertheless, we apply various kinds of resummation 
techniques [p4|, in order to obtain reliable quantitative results for the problem under consideration and to check the 
stability of these results. 

A. Chisholm-Borel resummation 

First, we employ a simple two- variable Chisholm-Borel resummation technique pl| . For our problem this turns out 
to be the most effective one. The resummation procedure consists of several steps: 

(i) starting from the initial RG function / in the form of a truncated series |34fl in the variables u and w, one 
constructs its Borel image: 

f-Va uV : y^ °ij("*) i M J " 

hi hJ 

where T(x) is the Euler's gamma function; 

(ii) the Borel image is extrapolated by a rational Chisholm (3^, |33| approximant [K/L] (ut, wt) which is defined 
as the ratio of two polynomials both in the variables u and w of degrees K and L such that the truncated Taylor 
expansion of the approximant is equal to that of the Borel image of the function /; 

(iii) the resummed function / rcs is then calculated as the inverse Borel transform of this approximant: 

00 
/ rcs = ( dtexp(-t)[K/L](ut,wt), 



There are a lot of possibilities to choose a Chisholm approximant in two variables. The most natural way is to construct 
it such that, if any of u or w is equal to zero, it leads to the familiar results for the reduced model. Here, for the 
Borel- images of the /3- functions Eqs. (|l|), ([19]) we have chosen the following approximants with linear denominators: 



/3 u (u,w,t) 



chis 



aifiut + CL2fiU 2 t 2 + aisuwt 2 + ao^w 2 t 2 + a2,iu 2 wt 3 + ai^uw 2 t 3 
1 + bifiUt + bo^wt 

\chis _ c 0,lUlt + C Q ^W 2 t 2 + C2,lU 2 wt 3 + C\^UW 2 t 3 



w (u,w,t) ct »* = ^ "" ' ;" —^ • (22) 

1 + difiut + do, lift 

Note, that the polynomials in the numerators are chosen to be symmetric in the variables u, w. In fig. we show 
the resummed 3d /^-functions in the t^w-plane for a = 2.9. In addition to the familiar fixed points describing 
Gaussian chains and polymers we obtain the stable LR fixed point for polymers in long-range-correlated disorder. 
For comparison, we depict the nonresummed functions in fig. || Only the Gaussian fixed point (u* = 0, w* — 0) is 
obtained without resummation. In figure ra we visualize the situation depicting the lines of zeroes of the resummed 
/3-functions at a — 2.9 in the u,w-plane in the region of interest. The intersections of these curves correspond to the 
fixed points. The corresponding values of the stable fixed point coordinates and the stability matrix eigenvalues for 
different values of the correlation parameter a < 3 are given in our Table [II] 

Substituting Eqs. (pp|), (fH|) into Eqs.(0)-(0) we get the following expressions: 



-if \ o l , (A ( a ) -/2(a)) . (c 2 {a) - Ci (a)) 23 2 n > 2 

v (u,w) = 2--uH uH UW + ——U + (c 3 (a) - c 1 (a))w , 

( , « 1 A(«) (/ 2 («)+4c 4 ( fl )-2A(a)) 1 2j _ 

(f 2 (a)-fi(a)f2(a)~Sc 3 (a)) 2 

16 [ ' 

This defines the critical exponents by v^ 1 = ^ _1 (w*, w*) and 7 = 7(1**, w*) at the stable accessible fixed point (u* , w*). 
To calculate these exponents in the region where the LR fixed point is stable, we again perform a resummation of the 
series in Eq. (|23|), using the following Chisholm approximants: 

cp.o + Cp^Wt + CifiUt + ci tl uwt 2 
1 + d\ out + do iwt 



The critical exponent r\ is obtained from the scaling law in Eq. (|12|). The numerical values for v, 7 and r\ are listed in 
Table In] for a — 2.3, . . . , 2.9. Note, that for a = 3, which corresponds to short-range-correlated point-like defects, the 
interactions u and w become of the same symmetry, so we pass to one coupling (u — w) and reproduce the well-known 
values of the critical exponents for the pure SAW model. The numerical values corresponding to those listed in table 
U are in this case: u* = 1.63, v = 0.59, 7 = 1.17, r\ = 0.02, uj — 0.64. As departing from the value a = 3 downward to 
2 one notices a major increase of the value of the coupling u, so the results are more reliable for a close to 3. At some 
value a — a marg the LR fixed point becomes unstable. This is explained by the following physical interpretation: as 
noted in the introduction, the case a = d — 1 (in our 3d approach a — 2) corresponds to straight lines of impurities 
of random orientation, and the absence of stable fixed points for a near a = 2 suggests the collapse of the polymer 
chain in such a medium. 

It is difficult to estimate the accuracy of the numerical values presented in Table M. On one hand, it is the first 
non-trivial result: in the one- loop approximation at fixed d, a one does not encounter the LR fixed point so one can 
not estimate deviations caused by different orders of the perturbation theory. On the other hand, as is known from 
the experience with the studies of magnets with long-range-correlated disorder [ pl| , the convergence of the resummed 
series for the RG functions is worse than in the pure (or the short-range-correlated) case. The resummed two-loop 
expansions we exploited here give quite reliable estimates for the exponents of SAWs on the pure (d = 3) lattice 
(compare our two-loop values v = 0.59 and 7 = 1.17 with the most precise RG estimates cited just after Eq. (|])) 



In the case of the short-range-correlated diluted magnets the comparison of recent six-loop results 1 35 with the two- 
loop ones [[3lj brings about the accuracy of latter of the order of several percents. While for our case no higher order 
calculations are available to test the numerical accuracy of the data in Table O the results clearly confirm the presence 
of a new stable fixed point LR with critical exponents that differ from those of the "pure" fixed point P. 

In order to confirm the quantitative stability of the picture we obtained, we have also used different non-symmetric 
approximants for /3 W instead of the one given in Eq. (|22|). As expected, this approach was less effective in the sense 
that the region where a stable fixed point could be established was reduced and the numerical values differed from 
those given in Table O. Nonetheless, the qualitative picture is the same: an LR fixed point exists and is stable in 
some interval a mar g < a < 3. 



B. Subsequent resummation 

Secondly, we applied the method of subsequent resummation, developed in the context of the d — dimensional 
diluted Ising model in ref. |3(| and successfully used for the d = 3 case in ref |3^]. Here, the summation was carried 
out first in the coupling u and subsequently in w. Starting from the /^-functions in Eqs. (Ilq) , (fill) we rewrite them 
as series in the variable w: 

2 95 3 / b% 2 



f3 u {u,w) = -u + u - — u +w I u(f 2 (a) - 3/i(a)) + — u 



:•■ / /■ , , ,.,- r^ , 6 e( a K \ , , ,, , ,,; 



+w 2 [ h(a) - «(6 8 (o) + -^) ) + & 5 (a)w 

We first perform a Pade-Borel resummation of the coefficients at different powers of w in the variable u, where it is 
possible. This results in a series of the form: 

f(u,w)=Y / M^W- (25) 



The coefficients Ai are some functions of u. Finally, the series (E5J) are resummed in the variable w. While we do not 
expect any high accuracy from this method, as far as the applicability has not been proven for our problem, again 
the presence of a stable fixed point LR for a marg < a < d in this case confirms the stability of a new type of critical 
behavior. 

We note that in addition to the above procedures we have tried a Pade-Borel approximation for the summation of 
the RG functions. To treat the two variable case we used the representation in terms of a resolvent series [[53| [36[ in 
a single auxiliary variable. However, no fixed points with i«/0 were found in the region of interest. But note, that 
even for the we akly diluted quenched Ising model this procedure does not lead to stable fixed points in the three loop 
approximation fbj]. 

C. Interpretation of the numerical results 



We may summarize and interpret our results as follows: 

(i) A new stable fixed point (LR) for polymers in long-range-correlated disorder is found for d — 3, a < d, leading 
to critical exponents that are different from those of the pure model; 

(ii) There is a marginal value a marg for the parameter a, below which the stable fixed point is absent, indicating a 
chain collapse of the polymer for disorder that is stronger correlated. 

(iii) The critical exponent v increases with decreasing parameter a, like in the Weinrib and Halperin case. But 
note, that the relation v = 2/a does not hold. Physically this means that in weak long-range-correlated disorder 
(a > a ma r g ) the polymer coil swells with increasing correlation of the disorder. The self avoiding path of the polymer 
has to take larger deviations to avoid the defects of the medium. 

V. CONCLUSIONS 

In the present work, we have analysed the scaling behavior of polymers in media with quenched defects that are 
correlated with a correlation that decays like ~ l/x a for large separations x. This type of disorder is known to 
be relevant in magnetic systems j20l fell , but the question about its relevance in the polymer problem was so far 
not answered. To this end we applied the field-theoretical RG approach, and performed renormalization for fixed 
mass and zero external momenta |25f. In our study we take special care of the symmetry properties of the effective 
Hamiltonian of the system H. Formerly |13| we performed calculations up to the linear approximation, using a double 
e, (5-expansion, as proposed for the magnetic problem in the work of Weinrib and Halperin J20J. While already this 
study indicated the possibility of a new type of critical behavior in such a system, it predicted such behavior for an 
unphysical range of parameters. A more sophisticated investigation at higher order of the perturbation series was 
needed to confirm the existence of a distinct polymer scaling behavior for long range correlated disorder. 



We use 2-loop expressions for the RG functions that were recently obtained for m-component systems in the fixed 
d, a approach |21j| and apply appropriate resummation techniques. This way we confirm that in a medium with 
long-range-correlated quenched disorder the swelling of the polymer coil is governed by a distinct exponent i^lr that 
increases when the correlation of the disorder is increased (i.e. a is decreased). When the correlation is too strong, 
i.e. a is below some marginal value a marg >~ 2, then a crossover to the collapse of the polymer is predicted. 

VI. APPENDIX 

Here, we turn our attention to the 3d magnetic system. Recently Ballesteros and Parisi J371] presented Monte- 
Carlo simulations of the site-diluted Ising model in three dimensions in the presence of quenched disorder with 
long-range-correlations. The values of the corrections-to-scaling exponents are of great interest in the interpretation 
of such simulations. In previous work, dedicated to 3d magnets with long-range-correlated disorder |16|, [17], |21j these 
exponents have not been calculated. Here we carry out these calculations based on the /^-functions of the model Eq. 
(0) in the 2-loop approximation, as presented in ref. pl[ . 

The correction-to-scaling exponent is defined as the minimal stability matrix eigenvalue in the stable accessible 
fixed point. We carry out the investigation in new variables (u, v, w) — ► (u, v,v + w), as proposed by Dorogovtsev M 



and perform the Pade-Borel resummation of the /3-functions. The results are presented in Table [II 



[1] J. des Cloizeaux, G. Jannink, Polymers in Solution. (Clarendon Press, Oxford, 1990); L. Schafer L, Universal Properties 

of Polymer Solutions as Explained by the Renormalization Group (Springer, Berlin, 1999). 
[2] P.-G. de Gennes, Scaling Concepts in Polymer Physics (Cornell University Press, Ithaca and London, 1979). 
[3] Recent renormalization group estimates for d — 3 critical exponents are found in: R. Guida, J. Zinn Justin. J. Phys. A 31, 

8104 (1998). 
[4] B. Nienhuis, Phys. Rev. Lett. 49, 1062 (1982); P. Grassberger, Z. Phys. B 48, 255 (1982). 
[5] B. K. Chakrabarti and J. Kertesz, Z. Phys. B - Condensed Matter 44, 221 (1981). 
[6] Y. Kim, J. Phys. C 16, 1345 (1983). 
[7] A.B. Harris, Z. Phys. B 49, 347 (1983). 

[8] Y. Meir and A. B. Harris, Phys. Rev. Lett. 63, 2819 (1989). 
[9] P. Grassberger, J. Phys. A 26, 1023 (1993). 
[10] K. Barat and B. K. Chakrabarti, Phys. Reports 258, 377 (1995). 

[11] S. B. Lee and H. Nakanishi, Phys. Rev. Lett. 61, 2022 (1988); S. B. Lee, H. Nakanishi and Y. Kim, Phys. Rev. B 39, 9561 

(1989); J. Machta and TR. Kirkpatrick, Phys. Rev. A 41, 5345 (1990); A.V. Izyumov and K. B. Samokhin, J. Phys. A 

32, 7843 (1999); A. R. Altenberger, I. I. Siepmann, and J. S. Dahler, J. Phys. A 272, 22 (1999); Y. Shiferaw and Y.Y. 

Goldschmidt, J. Phys. A 33, 4461 (2000). 

[12] For recent references see e.g. A. Ordemann, M. Porto, H. E. Roman, S. Havlin, and A. Bunde, Phys. Rev. E 61, 6858 

(2000); A. Ordemann, M. Porto, H. E. Roman, and S. Havlin, Phys. Rev. E 63, 0104 (2001). 
[13] V. Blavats'ka, C. von Ferber, and Yu. Holovatch, J. Mol. Liq. 91, 77 (2001). 
[14] About an influence of a non-correlated weak quenched disorder on a magnetic phase transition see e.g. R. Folk, Yu. 

Holovatch, and T. Yavorsk'ii, Phys. Rev. B 61, 15114 (2000) for a recent review. 
[15] A.B. Harris, J. Phys. C 7, 1671 (1974). 

[16] D. Boyanovsky and J. L. Cardy, Phys. Rev. B 26, 154 (1982). 
[17] S.N. Dorogovtsev, J. Phys. A 17, 677 (1984). 

[18] Y. Yamazaki, A. Holz, M. Ochiai, and Y. Fukuda, Physica A 150, 576 (1988). 
[19] Y. Yamazaki, A. Holz, M. Ochiai, and Y. Fukuda, Phys. Rev. B 33, 3460 (1985); Y. Yamazaki, Y. Fukuda, A. Holz, and 

M. Ochiai, Physica A 136, 303 (1986). 
[20] A. Weinrib and B.I. Halperin, Phys. Rev. B 27, 413 (1983). 

[21] V.V. Prudnikov, PV. Prudnikov, and A.A. Fedorenko, J. Phys. A 32, L399 (1999); V.V. Prudnikov, P.V. Prudnikov, and 
A.A. Fedorenko, J. Phys. A 32, 8587 (1999); V.V. Prudnikov, P.V. Prudnikov, and A.A. Fedorenko, Phys. Rev. B 62, 
8777 (2000). 
[22] G. Grinstein and A. Luther, Phys. Rev. B, 13, 1329 (1976). 
[23] E. Korutcheva and F. Javier de la Rubia, Phys. Rev. B 58, 5153 (1998). 

[24] E. Brezin, Le Guillou J C, Zinn- Justin J, "Field theoretical approach to critical phenomena" in Phase Transitions and 
Critical Phenomena, edited by C. Domb and M. S. Green, Vol. 6, Academic Press, London, 1976; D. J. Amit, Field Theory, 
the Renormalization Group, and Critical Phenomena (World Scientific, Singapore, 1989); J. Zinn- Justin, Quantum Field 
Theory and Critical Phenomena (Oxford University Press, 1996); H. Kleinert, V. Schulte-Frohlinde, Critical Properties of 
4> 4 - Theories (World Scientific, Singapore, 2001). 
[25] Parisi G., in Proceedings of the Cargese Summer School, 1973, (unpublished); J. Stat. Phys. 23, 49 (1980). 
[26] G. H. Hardy, Divergent Series (Oxford, 1948). 



10 

[27] J. P. Eckmann, J. Magnen, and R. Seneor, Commun. Math. Phys. 39, 251 (1975). 

[28] A. J. Bray, T. McCarthy, M. A. Moore, J. D. Reger, and A. P. Young, Phys. Rev. B 36, 2212 (1987). 

[29] A. J. McKane, Phys. Rev. B 49, 12 003 (1994). 

[30] G. Alvarez, V. Martin-Mayor, and J. J. Ruiz-Lorenzo, J. Phys. A 33, 841 (2000). 

[31] G. Jug, Phys. Rev. B, 27, 609 (1983). 

[32] J. S. R. Chisholm, Math. Comp. 27, 841 (1973). 

[33] G. A. Baker, Jn. and P. Graves- Morris, Pade Approximants (Addison-Wesley, Reading, MA, 1981). 

[34] We perform our analysis of the two-loop expressions for the RG functions in the fixed d, a technique. It is known that for 

other models with disorder, the £-expansion leads to less reliable results |l4J. 
[35] A. Pelissetto and E. Vicari, Phys. Rev. B 62, 6393 (2000). 
[36] P. J. S. Watson, J. Phys. A 7, L167 (1974). 
[37] H. G. Ballesteros and G. Parisi, Phys. Rev. B 60, 12912 (1999). 



11 



Fixed Point 


U* 


W* 




Wl 


W2 


Gaussian (G) 










— e 


-5 


Pure SAW (P) 


e 







£ 


e/2-<5 


Long-range (LR) 


2S A 
ie-S) 


S(s-2S) 


H*- 


- 4<J ± Ve 2 - 


-4e<5 + 8<5 2 } 



TABLE I: Fixed points and stability matrix eigenvalues in the first order of the e, 8 - expansion. 



a 


u* 


w* 


V 


7 


V 


Wl,2 


2.9 


4.13 


1.47 


0.64 


1.25 


0.04 


0.25 ± 0.62 i 


2.8 


4.73 


1.68 


0.64 


1.26 


0.04 


0.22 ± 0.76 i 


2.7 


5.31 


1.81 


0.65 


1.28 


0.03 


0.18 ± 0.89 i 


2.6 


5.89 


1.87 


0.66 


1.29 


0.03 


0.15 ± 0.99 i 


2.5 


6.48 


1.89 


0.66 


1.31 


0.02 


0.11 ± 1.09 i 


2.4 


7.10 


1.87 


0.67 


1.33 


0.01 


0.07 ± 1.18 i 


2.3 


7.76 


1.84 


0.68 


1.36 


0.01 


0.03 ± 1.26 i 



TABLE II: Stable fixed point of the 3d two-loop /3-functions, resummed by the Chisholm-Borel method, the corresponding 
critical exponents and the stability matrix eigenvalues at various values of a. 



a 


2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 


u){m = 1) 


0.80 0.81 0.83 0.87 0.94 1.14 1.07 0.87 0.71 0.69 


ui{m = 2) 


1.15 1.08 0.93 0.86 0.81 0.68 0.59 0.57 0.55 0.54 


ui{m = 3) 


0.88 0.83 0.76 0.67 0.62 0.61 0.60 0.60 0.59 0.68 



TABLE III: The correction-to-scaling exponents u> for the phase transition in the 3d m-vector model with m — 1, 2, 3. 




FIG. 1: The critical behavior of a polymer in a medium with long-range-correlated disorder in different regions of the d, a-planc 
as predicted by the first order e, <5-expansion. Region I corresponds to the Gaussian random walk behavior, in the region II 
scaling behavior is the same as in the medium without disorder, in region III the "long-range" fixed point LR is stable and 
the scaling laws for polymers are altered, in region IV no accessible stable fixed points appear; this may be interpreted as the 
collapse of the chain. 
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FIG. 2: The critical behavior of the m = 2-magnet in a medium with long-range-correlated disorder in different regions of 
the d, a-plane as predicted by the first order e, 5-expansion. Region I corresponds to the mean field behavior, in the region II 
the critical exponents are the same as in the medium with uncorrelated disorder, in III the fixed point LR is stable, in IV no 
accessible stable fixed points appear. 




FIG. 3: The Chisholm-Borel resummed 3d /3-functions in the two-loop approximation at a = 2.9. The flat surface corresponds 
to the /^-function. The resummation restores the presence of a pure SAW fixed point (u* — 1.63, w* = 0) and leads to a new 
stable "long-range-correlated" fixed point. The coordinate box is chosen to show the stable LR fixed point (u* = 4.13, w* = 1.47) 
on the face of the box. 
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FIG. 4: The non-resummed 3d /3-functions at a - 
the fixed points. Only the Gaussian fixed point u* 



2.9. The intersection of the U and f3 w surfaces with the u — w plane give 
= w* = is present without resummation. 




FIG. 5: The lines of zeroes of the 3d /3-functions (|l8[), ( |l9| ) resummed by the Chisholm-Borel method at a — 2.9. The dashed 
line corresponds to /3 U = 0, the solid lines depict (3 W — 0. The intersections of the dashed and solid lines give three fixed points 
shown by filled circles at u* — 0, w* — (G), u* = 1.63, w* — (P), and u* — 4.13, w* = 1.47 (LR). The fixed point LR is 
stable. 



